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ABSTRACT 


We numerically study the detailed evolutionary features of the wave-like disturbance and its prop- 
agation in the eruption. This work is a follow-up to Wang et al., using significantly upgraded new 
simulations. We focus on the contribution of the velocity vortices and the fast shock reflection and 
refraction in the solar corona to the formation of the EUV waves. Following the loss of equilibrium 
in the coronal magnetic structure, the flux rope exhibits rapid motions and invokes the fast-mode 
shock forward of the rope, which then produces the type II radio burst. The expansion of the fast 
shock, which is associated with outward motion, takes place in various directions, and the downward 
expansion shows the reflection and the refraction as a result of the non-uniform background plasma. 
The reflected component of the fast shock propagates upward and the refracted component propagates 
downward. As the refracted component reaches the boundary surface, a weak echo is excited. The 
Moreton wave is invoked as the fast shock touches the bottom boundary, so the Moreton wave lags the 
type II burst. A secondary echo occurs in the area where reflection of the fast shock encounters the 
slow-mode shock, and the nearby magnetic field lines are further distorted because of the interaction 
between the secondary echo and the velocity vortices. Our results indicate that the EUV wave may 
arise from various processes that are revealed in the new simulations. 

Subject headings: Sun: coronal mass ejections (CMEs) - Sun: magnetic fields - Magnetohydrodynamics 
(MHD) - Shock Waves 
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1. INTRODUCTION 


The most intense energetic activity in the solar sys- 
tem is the solar eruption that produces solar flare, erup- 
tive prominence, and coronal mass ejection (CME). Dur- 
ing the process, a large number of magnetized energetic 
plasmas (with mass of up to 1016 g and energy of 10°? 
erg) are ejected into the interplanetary space within a 
short timescale, and hence disturb spatial and planetary 
magnetic field and significantly affect satellite operation, 


aviation power, human space ex poration communica- 
tion and so on (Schwenn 2008): E (2007); [Lin] 
(2007); (2011); [Cheng et al} (2012) and references 


therein). One interesting “phenomenon closely associ- 
ated with CMEs is the globally propagating wave-like 
Gisturbanee | in the corona, i.e. the EIT or EUV wave. 
1| (1997) first reported this phenomenon, and 
(1998) first analyzed this phenomenon 
in detail using ‘the data from the Extreme-ultraviolet 
(EUV) Imaging Telescope (EIT) on (space) board the 
Solar and Heliosphere Observatory (SOHO) spacecraft. 

The EUV waves have generally been observed as broad, 
diffuse arc-shaped bright front with lifetime of about 50 
minutes. It could be seen in the lower corona (at temper- 
ature of 1-2 MK). The speed of the EUV wave front var- 
ied from about 50 to over 700 km s~ with ‘typical’ speed 
of 200-400 km s~! (Thompson & Myers|2009). me EUV 


waves are usually associated with EN dinm gs, type 


II radio bursts, and flares (Biesecker et al] 12009. Re- 
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Liu et al] (2011); [Zheng et al. 


cently, [Nitta et al] [Nitta ct al] (2013 (2013) presented a large sample of 
events that look like EUV waves as observed by the At- 
mospheric Imaging Assembly (AIA; [Lemen et al] (2019) 
on board the Solar a Observatory (SDO), and 
revisited their associations with flares, CMEs and type II 
radio bursts. They found that the speed of EUV waves is 
not strongly correlated with CME magnitude or the flare 
a. nor do they show an association with type II 
bursts (cf. [Nitta et al] (2014). 

a so far have shown that the EUV 
waves may largely be explained in terms of coronal 


2009): [Schmidt all Ofma: ‘Long et al! OU: 
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2 ; )) driven = ay: 
related eruptions although other components may 
co-exist in the same event or exist in other events 


(2012). It is possible to 

identify the properties of EUV waves that can be 

naturally linked with fast-mode shock waves. Obser- 

vations showed that in their early stages EUV mans 

either experience significant decelenat deceleration {heng et al. 
201% Bhen & Lin) 012a 201 


a a approximatel 


constant constant speede 


of a large number of samples. 

Regardless of their initial speeds or deceleration pro- 
files, these waves end up travelling within a speed range 
of 180-380 km s~!, which is consistent with the fast- 


2 Wang et al. 


mode speed over the quiet Sun 
[Patsourakos & Vonrlidas 2013). 

In addition, two other characteristics of EUV waves are 
found in case studies. First, what appears to be a single 
wave may in fact consist of fast and slow components. 
In the 2010 April 8 event, identified the 
two components and explained the large-scale ripples in 
terms of the faster component overtaking the slower one. 
Second, EUV waves are seen to reflect and refract as the 
propagate in the corona. See, for example, Shen et all 
(2013), who studied the 2012 April 23 event. Further- 
more, noticed in the 2011 August 4 
event that a secondary wave was excited by a reflected 
wave. 

Although the fast-mode MHD wave scenario of 
the EUV waves is supported by a lot of observa- 
tional and theoretical results (Thompson et_al. 
2004allb; 2009; 
Patsourakos & Vourlida 2009: 
SAF Loue ptala Bar Tiet 2011} [Li et_al] (201% 
Zheng et_al.)|2011) [2012 2012 
2012allb: , there is little understanding 
of how it evolves once the fast-mode shock is produced 
by solar eruptions. Therefore, numerical experiments 
are needed to better understand the details of how the 
EUV waves evolve. We demonstrate it here as our main 
goal of this work. 

In this paper, we numerically study the detailed wave- 
like disturbance and its propagation to uncover processes 
not commonly discussed but potentially important to un- 
derstand some aspects of the EUV waves. In our simula- 
tions, the employment of the high resolution grid and the 
ao eee d 
(1999), hereafter S&G) provides a good opportunity to 
study some important features that were not shown in 
our previous work 2012). 
We describe the physical model, formulae and numerical 
approaches in the next section. In Section 3, the numer- 
ical results are presented, and we discuss these results in 
Section 4. Finally, we summarize this work in Section 5. 
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2. NUMERICAL MODEL AND FORMULAE 


We suppose a two-dimensional magnetic configura- 
tion in the semi-infinite x-y plane. In the coordinates, 
y = 0 corresponds to the boundary between the photo- 
sphere and the chromosphere. The chromosphere and the 
corona are represented by y > 0. In this model, we use a 
current-carrying flux rope to represent a prominence or 
filament that floats in the corona, and a line dipole below 
the photosphere to denote the photospheric background 
field. The evolution of the magnetic system satisfies the 
following ideal MHD equations: 


where B indicates the magnetic field, J the electric cur- 
rent density, v the velocity of the flow, p the mass den- 
sity, p the gas pressure, e the internal energy density, G 
the gravitational constant, Mo the solar mass, Rọ the 
solar radius, y the ratio of specific heats, m, the pro- 
ton mass. The ZEUS-2D MHD code 
[1992allb; is employed to solve Equa- 
tions (1). 

In our simulations, the magnetic configuration is com- 
posed of the current-carrying flux rope, the image of the 
current inside the flux rope, and the background mag- 
netic field that is generated by a line dipole located at 
y = —d below y = 0. The relative strength of the 
dipole field M is denoted by a dimensionless parameter 
M = m/(Id), which is related to the ratio of the strength 
of the dipole field m and the product of the filament cur- 
rent J and the depth d of the dipole field 
[Wang et al] [2009). The flux rope is located at y = h 
above the bottom boundary. 

The initial magnetic configuration from which the 
eruption starts is: 


By = Bg(R_)(y — ho)/R- — Bo(R+)(y + ho)/R+ 

— Bg(r + A/2)Md(r + A/2)[a? — (y + d)?]/R4,(2) 
By =—Bg(R_)ax/R_ + Bg(Ry)x/Ry 

— Byg(r + A/2)Md(r + A/2)22(y + d)/ RG, (3) 


with 


R4 =y T (y = ho)”, 


R=" Tv (y T ay, 


and By(R) is determined by the electric current density 
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distribution j(R) inside the flux rope. They are: 


2 
B,(R) = —joR, 0<R<r—-A/2; 


sanel- S) (E) a de 


ORCE 


A 2 
+ (2) cos [E(R-r+3)|br-apenersare 
Bo(R) = -Zh [r? + (4/2)? — X(A/7)?] ,r + A/2 < R < o; 


2 


i(R)=0,r +4/2< R< œ. 


(R) 
j(R) = 2cos[a(R -r + A/2/A]+1,r - A/2 < R<r+A/2, 
j(R) 


(4) 


The initial background plasma density po(y) in this 
work is based on the S&G model, and therefore more re- 
alistic than that used in the previous model 
[Mei et _al][2012). We have: 


poy) = poof (y), 
Fly) = a12? (y)e"? [1 + as2(y) + aaz? (y) + a52°(u)] 5) 


z(y) = Ras 


Here we adopt poo = 1.672 x 10718 g cm~, which is 
about one order of magnitude smaller than that in our 
pre wus work. (Wans cad. and y denotes the 
height from the surface of the Sun. We take aj = 
0.001292, a2 = 4.8039, ag = 0.29696, a4 = —7.1743, 
ds = 12.321. 

The S&G model, which is two-dimensional and semi- 
empirical, was developed using the Skylab white-light 
coronagraph (Guhathakurta et alJ[1996) and Ulysses in- 
situ (Phillips et al. measurements. It smoothly 
connects the density distributions near and far from the 
Sun (see Figure 1 offLinl (2009). The density distribution 
f(y) given by Equation (5) describes an isothermal atmo- 
sphere. Density decreases exponentially with height in 
the lower corona, and then decreases much more slowly as 
y °. The results of radio observations of type II bursts 
over a wide frequency band from a few kHz to 13.8 MHz 
also suggested the y7? variation of the plasma density 
far Prout the oom (eeinne a ‘Tis moana (heh 
the atmosphere model with density distribution given in 
Equation (5) could be considered as a realistic model. In 
addition, we choose to use the present values of coeffi- 
cients for the density profile that describes the density 
distribution in the polar region. 

S&G also gives another profile for the density distribu- 
tion in the equatorial region where the helmet streamer is 
located (see equation [18a] and Table 1 of S&G’s original 
paper). Although the coefficients for these two profiles 
are different, the values of the resultant density distribu- 
tions in the two cases are not apparently different (see 
Figure 2 of S&G). This means that using either profile 
in our calculations would not affect the final results and 
conclusions in a major way. Figure [I] depicts the distri- 
butions of the initial configuration of the magnetic field 
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Figure 1. Distributions of the initial configuration of the mag- 


netic field (black contours) and the plasma density (shadings). The 
right color bar represents values of the density in lg p (g cm~*). 


(black contours) and the plasma density (shadings) used 
in the present work. 

We need to note here that the plasma density on the 
coronal base (y = 0.0) adopted in this work is about an 
order of magnitude higher than that in reality, but the 
magnetic field strength is roughly comparable to what 
is usually assumed. As a result, the Alfvén and related 
speeds in our work are about a factor of three lower than 
those that commonly appear in the literature. A lower 
Alfvén speed allows the fast shock to form more easily 
since we in this work aim to study the response of various 
layers of the solar atmosphere to the fast-mode shock 
driven by the lift-off of the flux rope. This may cause 
slower propagation and weaker strength of the shock, but 
will not change other properties of the shock. 

There is a balance between pressure gradient of gas 
and the gravity for the initial background atmosphere as 
below: 


GMo 
(Ro +y)? 
Taking Equations (5) into (6), we can get the initial back- 


ground pressure po(y). Then the temperature distribu- 
tion To(y) is obtained as follows 


Vpo(y) = —poly) (6) 


po(y) = LP kni), (7) 


where k is the Boltzmann constant. 

The initial total pressure consists of the gas pressure 
and the magnetic pressure. The initial total pressure and 
the mass density can be written as 


p= Bo(R)j(R)dR, (8) 
p= po(p/po)'/7. 


The computational domain is (—4L, 4L) x (0,8Z) with 
L = 10° km, and the grid points is 800 x 800. The bottom 
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Figure 2. The fast-mode speed cy contours at t = 0 s and t = 50s. 
The lower color bar represents values of the speed in km s~!. 


boundary at y = 0 applies a line-tied condition, while 
the other three use the open boundary. Table [I] lists the 
initial values of the parameters in our simulation. 


3. RESULTS OF NUMERICAL EXPERIMENTS 


In this section, we present the results of our numerical 
experiments. When the magnetic compression in our ex- 
periments surpasses the magnetic tension, the flux rope 
quickly rises from the beginning. Some important con- 
comitant phenomena appear as the flux rope rises, such 
as the fast-mode shock around the flux rope, the slow- 
mode shock emerging from two sides of the flux rope, 
and the velocity vortices. In order to let us more eas- 
ily find the region where the fast-mode shock occurs, we 
plot the distributions of the fast-mode magneto-acoustic 
speed, cf, at time t = 0 s and t = 50 s, respectively, as 
the eruption is in progress. We here have c+ = v4 + c? 


with va = |B|/\/47p and cs = \/yp/p being the Alfvén 
speed and sound speed, respectively. 

We studied disturbances caused by the interaction of 
the fast-mode shock with the boundary to explain the 
Moreton wave, and interaction between the slow shock 
and the velocity vortices to account for the EUV waves 
{Wang et all[2009). In this work, we focus on the forma- 
tion and propagation of the reflection/refraction of the 
fast-mode shock and the velocity vortices to investigate 
their contribution to the formation of the EUV waves. 

found that in the lower Alfvén speed 
environment, the velocity vortices and the slow shock 
may be responsible for EUV waves. In the higher Alfvén 
speed environment, the second echo of the fast-mode 
shock together with the slow shock and the velocity vor- 
tices may produce EUV waves. However, our results 
show that even in the same Alfvén speed environment, 
various origins of waves could be easily recognized due 
to the very high grid resolution used in our simulations. 


3.1. Formation of the vortices and the fast shock 


The loss of equilibrium results in fast upward move- 
ment of the flux rope and formation of a region with lower 
preseure bebiidiihe func top 
2009). This causes two flows. They are toward the lower 
pressure region, fetching in both magnetic fields of op- 
posite polarity and plasma. Forcing magnetic fields of 
opposite to flow together consequently leads to magnetic 
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Figure 3. The velocity divergence V -v and streamlines at t = 10 s 
and t = 20 s. The lower panels provide detailed structures in two 
subregions surrounded by rectangles marked in the upper panels, 
respectively. The blue cross shows the location of the left edge of 
the flux rope, and the red cross denotes the location of the velocity 
vortices at 10 s. The unit of time is second. The right color bar 
represents values of the velocity divergence in arbitrary unit. 


reconnection in this region. The velocity vortices on ei- 
ther side of the reconnection region forms as shown in 
Fig.[3] 

To look into more details of the velocity vortices and 
their impact on the flux rope, we enlarged the region of 
interest and created a composite of the velocity diver- 
gence V - v and streamlines in Fig.B] In the lower left 
panel, the blue and red crosses show, respectively, the lo- 
cations of the left edge of the flux rope and the velocity 
vortices at 10 s, respectively. Figure [4| shows the loca- 
tions of the velocity vortices on the left side of the flux 
rope (in red) and the left edge of the flux rope (in blue) 
at various times. 

From Figs.Bland[4] we notice that the velocity vortices 
blend with the edge of the flux rope at t = 10 s. It is not 
easy to make a distinction between the velocity vortices 
and the edge of the flux rope. However, they start to 
separate, and the distance between them increases with 
time. 

Another important finding is that a fast-mode shock 
starts to form in front of the flux rope at the time when 
the velocity vortices almost separate from the flux rope. 
To reveal more details of this process, we check the for- 
mation time of the fast-mode shock as shown in Fig.[5] 
It plots variations of the plasma density versus height 
along the y-axis at two different times. The highest peaks 
in both panels indicate the flux rope. In the panel of 
t = 10 s, we can recognize a slight increase in the density 
right forward of the flux rope, which represents the for- 
mation of the fast-mode shock. The fast shock is more 
clearly seen in the panel of t = 20 s, which may corre- 
spond to the occurrence of type II radio bursts as dis- 
cussed by [Wang et al] and [Lin et al] (2006). Just 
during this period, the velocity vortices begin to separate 
from the edge of the rope. 
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Table 1 
Initial values for several important parameters of the numerical experiment 


ro = 2.5 x 10° km ho = 6.25 x 10° km 


poo = 1.672 x 10713 g cm~? 


red: location of vortices 
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Figure 4. The location of the velocity vortices on the left side of 
the flux rope, and the location of the left edge of the flux rope at 
various times. 
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Figure 5. Plasma density distribution along the y-axis at t = 10 s 
and t = 20 s. The highest peak in two panels represents the flux 
rope. In the panel of t = 10 s, a slight increase indicates the 
formation of the fast-mode shock, and in the panel of t = 20 s, the 
fast shock is formed. 


Figure[6]plots the evolutions of V -v and streamlines as 
the flux rope moves outward. Three regions surrounded 
by rectangles display places where the flow is strong. De- 
tails in these three regions can be seen more clearly in the 
lower enlarged panels in Fig.[6] The plasma flow man- 
ifests apparent vortices near the side back of the flux 
rope. In each panel, a crescent feature around the flux 
rope represents the fast-mode shock. With the fast shock 
propagating forward, it expands sideward and backward, 
and its footprints touch the regions of the vortices at 
about t = 20 s. To obtain more details of this process, 
we studied the behavior of the distribution of V -v on 
the layer of the vortices, i.e. y = 1.4 x 104 km. Ten 
curves in Fig.[fJare for the distributions of V - v on this 
layer (y = 1.4 x 10+ km) at different times. Plots in 
Fig.[7] clearly show the fast shock that sweeps the layer 
y = 1.4 x 104 km from t = 20 s, but there is no sign of 
the fast-mode shock before t = 10 s. The corresponding 
speed of the fast shock in this layer is about 300 km s71. 
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Figure 6. Evolutions of the velocity divergence V -v and stream- 
lines at various times. The lower panels provide detailed structures 
in three subregions surrounded by rectangles marked in the upper 
panels, respectively. The right color bar represents values of the 
velocity divergence in arbitrary unit. 


Figure [3] suggests that the footprint of the fast shock 
does not arrive at the layer where the vortices occur un- 
tilt = 20 s. Our simulations reveal that the impact of the 
fast shock joins that of the velocity vortices from t = 20 s, 
inducing more complex patterns of flow or disturbance in 
this layer. On the other hand, the fast-mode shock it- 
self does not stay in this layer. Instead, it continues to 
expand outward, leaving behind the disturbances that it 
causes in the corona. Obviously, the fast shock propa- 
gates in the corona faster than the disturbances that it 
produces. This might account for the fast component of 
EUV waves usually observed in the same coronal layer. 


3.2. Formation and propagation of the fast shock 
reflection and refraction 


In addition to the velocity vortices and the fast-mode 
shock, one more feature is shown in Figs.[6] and B] A 
composite of the plasma density (shadings) and magnetic 
field lines (black contours) during the interval between 
50 s and 260 s can be seen in Fig [8] At t=50 s, the 
plasma flow near the side back of the flux rope results in 
vortices, as clearly seen in the lower-left panel of Fig.[6] 
The flow is so strong that it impacts the nearby magnetic 
field lines, resulting in distinct deformation of the nearby 
magnetic structures. Because the background plasma in 
the corona is not uniform, reflection and refraction take 
place at both wake ends of the fast shock as it moves 
downward (see Fig.[8). The reflected component of the 
fast shock propagates upward and the refracted compo- 
nent propagates downward continuously. In this paper, 
we use the term “reflection” for the reflection that takes 
place in the propagation path of the fast-mode shock due 
to the non-uniform media. The “echo” for the reflection 
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Figure 7. Evolution of V -v on the layer y = 1.4 x 10* km at 
various times. 


occurs on the boundary surface. However, in physics, 
reflection and echo are actually the same phenomenon. 

With the evolution progressing, two wake ends of the 
fast-mode shock sweep side-back to the flank of the flux 
rope as shown in Fig.[8] The incident angle of the fast 
shock at the layer of its wake ends becomes large. Pan- 
els in Fig.[8] clearly show that the angle of the reflection 
and the angle of the refraction vary with the incident 
angle. To further study this point, we investigate the 
variation of the angle of the reflection and the angle of 
the refraction as the flux rope moves outward. Figure [9] 
plots the evolutions of the angles of incidence, reflection 
and refraction, respectively, at t = 70 s and t = 190 s. 
With the downward-propagating fast shock encountering 
denser background plasma, the incident angle becomes 
large, which results in the increasing reflection and re- 
fraction angles between t = 70 s and t = 190 s. 

We also investigate the refractivity of the fast shock 
at its wake ends and the distribution of the background 
plasma density versus height as shown in Fig.[10] Solid 
points indicate the fast shock refractivity at different 
times. The curve represents the background plasma den- 
sity distribution versus height, normalized at y=0. From 
Fig. [10] we can see that the background plasma density 
and the refractivity increase with decreasing height. 

From the panels in Fig.[8] we see that the refraction 
of the fast-mode shock propagates downward as the fast 
shock extends outward and sideward. At about t = 100s, 
the refracted wave touches the boundary surface, and 
then reflection of the refracted wave occurs at the bound- 
ary surface. An echo of the refracted wave is thus pro- 
duced at each of the refracted wave’s wake ends, which 
goes back into the corona. This echo is weak, so it can 
not be easily distinguished. At about t = 260 s, the fast- 
mode shock sweeps the boundary surface. Because the 
background density gradient becomes sharp, the strong 
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Figure 8. Set of snapshots of magnetic field (black contours) and 
the plasma density (shadings) interval between 50 s and 260 s. The 
right color bar represents values of the density in lgp (g cm™?). 
Two regions surrounded by rectangles at t = 70 s and t = 190 s 
display areas where the angle of the reflection and the angle of the 
refraction vary with the incident angle growing large. Details in 
these two regions can be seen more clearly in Figure P] 
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Figure 9. Evolutions of the angle of incident, the angle of reflec- 
tion and the angle of refraction at two different time. a1 and a2 
correspond to the angle of incident at t = 70 s and t = 190 s, bı 
and bə are for the angle of reflection at t = 70 s and t = 190 s, cy 
and c2 denote the angle of refraction at t = 70 s and t = 190 s. 


reflection of the fast shock appears. 

Figure plots 29 distributions of V - v curves ver- 
sus x on the boundary y = 0 in the time interval be- 
tween t = 20 s and t = 580 s, which gives more in- 
formation about propagation of the refracted wave and 
the fast-mode shock at the boundary surface. The plots 
clearly display two different features of propagation at 
the boundary: the red arrow is for the footprint of the 
refracted wave and the blue arrow is for the footprint 
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Figure 10. The distribution of the background plasma density 
f(y) and variation of the fast shock refractivity versus height. Solid 
points indicate the fast shock refractivity at different time. The 
curve represents f(y) versus height. 
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Figure 11. Distributions of V -v on the boundary layer y = 0.0 
in the time interval between 20 s and 580 s. 


of the fast-mode shock. Figure [LI] further confirms that 
the refracted wave approaches the boundary surface at 
about 100 s. The fast shock starts to sweep the surface 
at about 260 s as manifested in Fig.[8} By measuring the 
distances that they propagate in the given time interval, 
we find that the speeds of the refracted wave and fast 
shock are 93 km s~! and 165 km s71}, respectively, when 
they sweep the boundary surface at y = 0. 

As demonstrated in Figs.[6] and [8] the fast shock, its 
reflected wave and refracted wave (or reflection of the 
refracted wave) interact with the vortices. At this stage, 
their interaction could be related to the origin of the 
EUV waves. 


3.3. Formation of the Moreton wave and the secondary 
echo of the fast shock 


With the fast shock propagating outward and expand- 
ing downward, the Moreton wave is invoked as the fast 
shock touches down at about t = 260 s. As we illus- 
trated in subsection 3.1, the fast shock forms in front of 
the flux rope at about t = 20 s, which then produces the 
type TI burst [Lin et al]|2000). This 
result indicates that the Moreton wave lags the type H 
burst about 4 minutes, which is helpful to understand 
their correlation observed in solar eruptions. 

As the fast shock touches the boundary surface, the 
echo of the fast shock will occur because the background 
density gradient becomes sharp. The echo at the bound- 
ary mixes with the reflection of the pre-existing wave, 
so the echo that propagates back into the corona be- 
comes so strong. It moves sideward much faster than 
the initial reflected wave (i.e. it appears on the layer 
above the boundary surface, rather than the surface) of 
the fast shock. Meanwhile, the initial reflected wave de- 
cays significantly since the footpoint of the fast shock 
moves outward and downward. On the other hand, the 
velocity vortices propagate sideways and upward, so the 
impact of the initial reflection of the fast-mode shock on 
the vortices becomes weaker and weaker as illustrated in 
Fig.[2] In this figure, we plot the evolutions of V -v and 
streamlines at 300 s and 400 s. The lower panels provide 
detailed structures in two subregions surrounded by the 
rectangles as marked in the upper panels. 

The initial reflected wave of the fast-mode shock decays 
gradually and becomes weak. However, the echo of the 
shock is fast and strong, producing one large structure 
with two small vortices when it enters the region where 
the vortices and the initial reflected wave are located. 
One of the two small vortices is influenced by the initial 
reflected wave that appears at the same layer above the 
boundary surface, while another one is impacted by the 
echo of the fast shock. The latter is stronger and faster 
than the former, overtaking the former and producing 
multiple “ripples”. observed similar 
phenomena in the 2010 April 8 event. 

With the flux rope moving up, the slow-mode shock 
is also generated in front of the flux rope as shown in 
Fig.13} This figure displays evolutions of V -v and 
streamlines at various times. Unlike the fast-mode shock, 
the slow shock cannot reach the boundary surface, be- 
cause it starts to decay and dissipate above the surface 
. From Figs.[3] and 
[14| we see that the echo of the fast-mode shock comes 
across the slow-mode shock before the slow shock dis- 
sipates, and the reflection of the fast shock’s echo (i.e. 
the secondary echo of the fast shock)is excited in the re- 
gion where the echo of the fast shock encounters the slow 
shock as shown in Fig.{I4] 

We also investigate the distributions of the z- 
component of V xv at layer of 3x 10+ km from the surface 
of the Sun as we did in the previous work 
(2009). Figure plots such distributions at various 
times. The curve at the bottom is for t = 50 s, and 
the one at the top is for t = 800 s. From Fig.[I5] we can 
see four propagating signs: the blue arrow corresponds 
to the fast shock; the red arrow is for the refracted wave 
of the fast shock; the green arrow shows the echo; the 
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Figure 12. Evolutions of V -v and streamlines at t = 300 s and 
t = 400 s. The lower panels provide detailed structures in two 
subregions surrounded by rectangles marked in the upper panels, 
respectively. The right color bar represents values of the velocity 
divergence in arbitrary unit. 


black one denotes the secondary echo. 

As we mentioned before, the inhomogeneity of the 
plasma causes the fast shock to reflect and refract in the 
way as the shock propagates towards the bottom bound- 
ary. The shock reaches the layer of 3 x 104 km at about 
t = 100s. The reflection and the refraction of the shock 
take place simultaneously at the same layer. At about 
t = 250 s, the fast shock arrives at the bottom boundary 
where the reflection becomes the echo and the refracted 
wave disappears in our calculation. In fact, as indicated 
in Fig.[8] the refracted wave should arrive at the bottom 
boundary earlier than the fast-mode shock itself and pro- 
duces an echo as well. But this echo is weak and blended 
with other kinds of disturbances in the corona, so we 
cannot recognize it in the present simulation. 

After t = 250 s, the fast shock sweeps the bound- 
ary surface, and its echo propagates back to the corona. 
At about t = 300 s, the echo gets to the layer of 
y = 3 x 10+ km, which serves as an extra source of dis- 
turbance in this layer as indicated in Fig.{[5} The echo 
meets the slow mode shock at t = 500 s, which gives rise 
to a further disturbance in the region. Consistent with 
our previous discussions about Fig.[7] the disturbance di- 
rectly caused by the fast shock itself always moves faster 
than the others that are caused by various types of sec- 
ondary effect and their combinations. 

By measuring the features seen in Fig.[15} as we did in 
Fig [Td] we find that the speed of the fast shock at layer 
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Figure 13. Evolutions of V -v and streamlines at various times. 
The lower panels provide detailed structures in three subregions 
surrounded by rectangles marked in the upper panels, respectively. 
The right color bar represents values of the velocity divergence in 
arbitrary unit. 


y = 3 x 10* km is 202 km s7}, which is slightly slower 
than that at a lower layer (cf. Fig.[Z). In addition, the 
speeds of the refracted wave, the echo and the second 
wave, at the same layer, are 147 km s71, 123 km s7}, 
and 50 km s~!, respectively. Furthermore, Fig.[15] de- 
picts a scenario of wave propagation in the corona where 
the EUV waves are usually observed such that a fast 
component is followed by a group of slower components 
with a large variety. 

The interaction between the secondary echo propagat- 
ing downward with the velocity vortices results in dis- 
tortion of the magnetic field line as indicated in Fig.|16} 
which shows the evolutions of magnetic field and plasma 
density. The continuous contours show magnetic field 
lines and the color shading the density distribution. Fig- 
ures and [16] indicate that the EUV waves in the 
last stage of propagation may be contributed by combina- 
tions of the echo and the secondary echo of the fast-mode 
shock, the velocity vortices, and the slow shock. 


4. DISCUSSIONS 
Following our previous work (Wang et al] , we 


further studied detailed evolutionary features of wave- 
like disturbances and possible relation and/or contribu- 
tions of these features to the formation of the EUV waves. 
In our study, the S&G empirical atmosphere model is 
employed for the density distribution of the background 
field, and the present simulation has much higher reso- 
lution than that of (2009). With these im- 
provement, we are able to perform more investigations on 
the object of interest and to inspect more aspects of the 
relevant issues. However, we need to note here that the 
purpose of our present work is not to reproduce a specific 
EUV wave event as the works done b 
and (Downs et al] (201 2012). Instead, our em- 
phasis is to investigate how the solar atmosphere would 
respond to an eruption that produces CME and other 
associated activities, to explore what observational con- 
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Figure 14. Evolutions of V -v and streamlines at t = 1200 s. 
The echo, secondary echo, slow shock, and vortices are denoted in 
the upper panel. The lower panel shows detailed structure in a 
subregion surrounded by rectangle marked in the upper panel. 


sequences we could expect in such processes, and to see 
how these consequences would contribute to the forma- 
tion of EUV waves. 

The most important discoveries or results of this work 
are the reflection and refraction of the fast-mode shock 
during its propagation, and formation of the secondary 
echo when the first echo of the fast shock interacts with 
the slow mode shock (Wang et_al][2009). No numerical 
experiments have ever reported such a phenomenon be- 
fore. The various disturbances that we identify in this 
work may correspond to actual observational features re- 
lated to the EUV wave phenomenon observed in many 
solar eruptive processes. Below we discuss our results 
one by one. 

As the flux rope goes up, a fast-mode shock commences 
to form forward of the flux rope at about t = 20 s, which 
then produced the type II radio bursts 
[2006). Another important finding is that the 
velocity vortices begin to separate from the flux rope 
when the fast shock forms. These results duplicate what 
observed in a specific event. 

Associated with the outward motion, expansion of the 
fast shock takes place in various directions. The down- 
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Figure 15. Distributions of (Vxv)z at the layer of y = 3x104 km 
that is in the corona. The time intervals for drawing these curves 
are between 50 s and 800 s. 


ward expansion shows the reflection and refraction as a 
result of the non-uniform background plasma density. In 
the uniform atmosphere (Forbed/990 [Wang et LEUN, 
the fast-mode shock propagates through it without reflec- 
tion and refraction. In the non-uniform atmosphere, on 
the other hand, reflection and refraction inevitably occur 
when the fast shock reaches the boundary, separating the 


media of different densities. Because the density changes 
continuously in the real corona, reflection and refraction 


take place everywhere in the propagation path of the 
shock. In our previous work (Wang et_al.|2009; 


2012), we did not notice this phenomenon because of the 


low resolution of the grid used in those calculations. 

The reflected component of the fast shock propagates 
upward, whereas the refracted one propagates downward. 
As the refracted component reaches the bottom bound- 
ary, an echo inevitably occurs. But this echo may not be 
strong enough, so we cannot distinguish its effect from 
that of other types of disturbance in the same region 
(see Figs.[6]and[8). However, it surely contributes to the 
formation of the slow component of EUV waves because 
it is always located behind the fast-mode shock itself. 
So our simulations show that the origin of EUV waves 
includes not only the fast shock propagating freely in a 
certain layer, but also various disturbances caused by the 
fast shock and the motion of the flux rope. This implies a 
complex origin of EUV waves. This scenario is consistent 
with the results of (2013). 

With the fast shock expanding backward sequentially, 
the background plasma becomes more dense. As a result, 
the incident angle becomes large, increasing the reflection 
and refraction angles. We also studied the refractivity 
of the fast shock versus height, which shows that the 
background plasma density and the refractivity increase 
with decreasing height. At about t = 100 s, the refracted 
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Figure 16. Evolutions of magnetic field (black contours) and the 
plasma density (shadings). The unit of time is second. The right 
color bar represents values of the density in lg p (g cm~*). A region 
surrounded by rectangle displays an area where the interaction be- 
tween the secondary echo propagating downward with the velocity 
vortices results in distortion of the magnetic filed line. Details in 
this region can be seen more clearly in the lower panel. 


wave touches the boundary surface and gets reflected, so 
an echo of the refracted wave is created at each of the 
refracted wave’s wake ends. This echo is weak, so it is 
not easy to be distinguished. 

When the fast shock sweeps the boundary surface, the 
strong echo of the fast shock is produced because of the 
sharp density gradient on the boundary. And then More- 
ton wave is invoked. The Moreton wave lags the type II 
burst that is produced by the formation of the fast shock 
about four minutes before, which is helpful to under- 
stand their correlation observed in solar eruptions. Us- 
ing the observations of the Hiraiso Radio Spectrograph in 
the frequency range of 25-2500 MHz, 
showed the observed type II radio burst is connected to 
the Moreton wave and fast-bright EUV wave. Further- 
more, their results also support that the Moreton wave 
originates from the fast shock. 


More details of these features are revealed by the ve- 
locity divergence V - v curves versus x on y = 0, which 
clearly displays two different features of propagation at 
the boundary surface: one is the footprint of the refracted 
wave; another is the footprint of the fast shock (Fig [II). 
The result shows that the refracted wave speed at the 
boundary is usually 1/2 that of the fast shock at the 
boundary. 

In addition to the fast-mode shock, the slow-mode 
shock can also be recognized, which is caused by the 
movement of the flux rope. Unlike the fast shock, the 
slow shock begins to decay above the boundary sur- 
face, so it cannot reach the surface. Because the echo 
of the fast shock from the boundary surface is strong, 
the secondary echo is excited when the echo of the fast 
shock encounters the slow shock at some layer above the 
boundary surface. This phenomenon may have been ob- 
served in the events reported by and 
(2013). 

Furthermore, we point out that discussions we made 
above revealed rich information about the origin of 
EUV waves. We noticed that, except the distur- 
bance caused directly by the fast-mode shock itself, 
all the other types of disturbance lag behind the 


fast shock, which has been confirmed by many ob- 
servations (e.g. ; [Asai et_al| (2012): 
(201 2allb)); (2013); 
Chen et al. explained the EUV 
waves behind the fast-mode shock as a result of the field 
line stretching by the CME motion, and therefore the 
EUV wave is not a true “wave”, but a so-called pseudo 
wave. Obviously, our results here, together with those 
of Wang et al] (2009) and Mai et all (2019), provide an 
alternative explanation for the results by 


(2005). 


5. CONCLUSIONS 


In this work, we use empirical atmosphere S&G model 
and a high resolution grid for numerical simulations to 
study detailed evolutionary features of wave-like distur- 
bance. Our main conclusions are summarized below. 

1. Following the loss of equilibrium in the coronal mag- 
netic structure, the flux rope exhibits rapid motions, and 
invokes the fast-mode shock in front of the rope, which 
then produces the type II radio bursts. 

2. The velocity vortices form on either side of the flux 
rope, and the vortices mix with the edge of the flux rope 
in the beginning. However, the vortices start to separate 
from the edge of the flux rope when the fast shock forms. 

3. The fast shock expands downward and sideward si- 
multaneously as it propagates forward. The downward 
expansion produces reflection and refraction as a result 
of the non-uniform background plasma density. The re- 
flected component of the fast shock propagates upward 
and the refracted component propagates downward. 

4. As the refracted component touches the boundary 
surface, a weak echo is produced. The Moreton wave is 
invoked as the fast shock sweeps the bottom boundary, 
so the Moreton wave lags the type II burst. 

5. As the fast shock touches the boundary surface, 
the strong echo occurs because the background density 
gradient becomes sharp. This echo propagates back into 
the corona. 

6. Besides the fast-mode shock, the slow-mode shock 
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is also excited by the movement of the flux rope. Un- 
like the fast shock, the slow shock begins to decay above 
the bottom boundary, so it cannot reach the boundary. 
The secondary echo occurs in the area where the echo 
of the fast shock encounters the slow shock. The nearby 
magnetic field lines are further distorted because of the 
interaction between the secondary echo and the vortices. 

7. During the subsequent evolution, we found an inter- 
esting phenomenon: a large structure consisting of two 
small vortices. One vortex is influenced by the initial 
reflected component of the fast shock, and another one 
is impacted by the strong echo of the fast shock. The 
latter vortex is stronger and faster, and so it overtakes 
the former in a short time scale, which may correspond 


to the ripples as observed in a specific event (e.g. see 
(2010). 


Our results indicate that various origins of the EUV 
waves exist. These include the fast shock, the reflected 
and refracted component of the fast shock, the echo of 
the refracted component, the echo and the secondary 
echo of the fast shock, the slow shock, and the ve- 
locity vortices (e.g., (1990); [Wang et al] (2009): 
Pen Considering that the EUV waves ap- 
pear roughly in the same layer as the disturbance that 
is usually left behind the fast-mode shock, we conclude 
that the slow component of the EUV waves may be in- 
voked by one of (or the combination of) these sources, 
and the fast EUV wave component should be due to the 
fast-mode shock itself. 
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